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ABSTRACT 

The large majority of extragalactic star cluster studies performed to date essentially use multi- 
colour photometry, combined with theoretical stellar synthesis models, to derive ages, masses, 
extinction estimates, and metallicities. M31 offers a unique laboratory for studies of globular 
cluster (GC) systems. In this paper, we obtain new age estimates for 91 M31 globular clus- 
ters, based on improved photometric data, updated theoretical stellar synthesis models and 
sophisticated new fitting methods. In particular, we used photometric measurements from the 
Two Micron All Sky Survey (2MASS), which, in combination with optical photometry, can 
partially break the well-known age-metallicity degeneracy operating at ages in excess of a few 
Gyr. We show robustly that previous age determinations based on photometric data were af- 
fected significantly by this age-metallicity degeneracy. Except for one cluster, the ages of our 
other sample GCs are all older than 1 Gyr. Their age distribution shows populations of young 
and intermediate-age GCs, peaking at ~ 3 and 8 Gyr respectively, as well as the "usual" com- 
plement of well-known old GCs, i.e., GCs of similar age as the majority of the Galactic GCs. 
Our results also show that although there is significant scatter in metallicity at any age, there 
is a noticeable lack of young metal-poor and old metal-rich GCs, which might be indicative 
of an underlying age-metallicity relationship among the M3 1 GC population. 
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1 INTRODUCTION 

M 31 is the nearest lar ge spiral galaxy, at a distance of about 780 
pc IStanek & Garnavichll998tlMacril200ll) . It contains more than 
337 confirmed globular clusters (GCs) and about 688 GC candi- 
dates (Galleti et al. 2004), i.e., significantly more than in our own 
Galaxy. Thus, M31 provides an excellent opportunity to study the 
properties of a large sample of GCs. From the observational ev- 
idence collected thus far (see, e.g.. iRich et al. I l2005t) . the M31 
GCs and their Galactic counterparts reveal some strikin g similari- 
ties iFusi Pecci et al. ll994llDioreovski et al. ll997tlBarmbv et al. I 
2002). For example, both GC systems seem to have similar mass- 
to-light ratios, structural parame ters, and ve locity dispersion - lu- 
minosity relations (see also lde Griis et al. l2005bl) . Studies of GCs 
in M3 1 can not only throw light on lots of questions about the for- 
mation, evolution and properties of M31 itself, including its mass, 
dynamics and chemical composition, they can also improve our un- 
derstanding jjfthe_formation and structure of galaxies in general 
(Battis tini et al. Il 980). In addition, GCs can provide us with good 
samples of Population II stars characterised by homogeneous abun- 
dances and histories, and with unique stellar dynamical conditions 
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for our study I Barmb v & Huchra 120011) . Therefore, GCs are so im- 
portant that they are considered as the fo ssils of the earliest stages 
of galaxy formation 1 Barmbv et al. 2000). However, at the distance 
of M31, construction of colour-magnitude diagrams (CMDs) be- 
low the main sequence turn-off, the most reliable method for age 
determinations of stellar populations, is extremely challenging for 
current state-of-the-art instrumentation. Here, one suffers from the 
dual effects of crowdi ng and the intrinsic faintness of the c luster 
stars jRich et al. Il996ll2005l:lstephens et al. 11200 lklBeaslev et al. I 
|2004 . although we note that lBrown et al. U2004I) presented a CMD 
down to the turn-off for a ~ 10 Gyr-old M31 GC. 



Since the pio neering work of iTinslev I 1 19681 1 19721) and 
ISearle et al. I ll973l) . evolutionary population synthesis modeling 
has become a powerful tool to interpret integrated spectrophoto- 
metric observations of galaxies as well as their components (see, 
e.g.. lAnders et al. 120041) . Comprehensive compilations of relevant 
curren t mode l sets, such as e.g. developed by Bruzual & Chariot 
ll993l Il99d henceforth BC93, BC96). Leitherer & Heckman I 

11997J), were provided by 
, The evolution of star 



1 1995), and Fioc & Rocca-Volmerange 
Leitherer et al. Ul996l) and lKennicuttNl998 
clusters is usually modeled based on the "simple stellar population" 
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(SSP) approximation 1 , which is a highly robust approximation for 
old G Cs in par t icular (e.g.lBarmbv & H uchra 20 0(J). 

iMa et al. I J200ll 12002 Jblcl) and l.Tiang et all <2003l) estimated 
the ages of, respectively, 180 star clusters in M33 and 172 GC 
candidates in M31 by comparing the SSP synthesis models of 
BC96 with the clusters' integrated photometric measurements in 
the Beijing-Arizona-Taiwan-Connecticut (BATC) photometric sys- 
tem. |MaetalJ y 006a) also determined the ages and metallicities 
of 33 M31 GCs and candidates using the updated method and up- 
dated SSP synth esis models of iBruzual & Chariot I 120031 hence- 
forth BC03), and Ma etaTl l2006bl) estimated the age and redden- 
ing value of the M31 GC 037-B327 based on photometric measure- 
ments in a large number of broad-band passbands from the optical 
to the near-infrared. 

From an observational point of view, the study of M31 
GCs is complicated, since in most cases we only have access 
to their integrated spectra and photometry, and cannot study 
the resolved stellar population. Therefore, we can only obtain 
the key physical parameters, such as the age and metallicity, 
by analysis of the integrated spectra or photometry. How- 
ever, a large body of evidence suggests that there is a strong 
age-metallicity degen eracy if on l y optical pho tometry is used 
lArimotol Il996t IWorfhevI 1 1994 iKavirai et all l2006h . A very 
useful method to break this degeneracy is through the application 
of particular spectral diagnostics based on the occurrence of 
individual st e llar absorption-line features (e .g.,|F^b^rll9_73|;|RosgJ 
1984 Il985t iDiaz. Terleyich & Terlevich I Il989t IWorthev et al I 
1994 IjonesfeWorthevI Il995t IVazdekis & Arimoto I Il999t 
Bruzual & Chariot 120031) ' 

At the same time, observational GC spectral energy distribu- 
tions (SEDs) are affected by reddening, an effect that is also dif- 
ficult to separate from the combined effects of age and me tallicity 
lCdzettilfl99"7l IVazdekis et"aLlll997t loriglia et al. 1119991) . How- 
ever, if the metallicity and reddening are derived accurately (and, 
ideally, independently), these degeneracies are largely (if not en- 
tirely) reduced, and ages can then also be estimated accurately 
based on a comparison of multi-colour photometry spanning a 
significant wavelength range Griis ct al. 2003b; Ander s et al. I 
2004) with theoretical stellar population synthesis models. 

In this paper, we present new age estimates for 91 GCs of 
the Ijiang et al. I 120031) sample, based on their 13 intermediate- 
band photometry in the BATC system, combined with additional 
broad-band optical and 2MASS near-infrared photometry, and on 
the updated SSP synthesis models of BC03. Section [2] describes 
the intermediate-band, broad-band and 2MASS photometry of our 
GC sample. In Section|3| we describe the metallicity and reddening 
data of the sample GCs; Section|4]includes a description of the SSP 
models used, and of our method to estimate the ages of the sample 
GCs. The main results and a discussion are also presented in this 
section. We summarise and conclude the paper in Section|5| 



2 INTERMEDIATE-BAND, BROAD-BAND AND 2MASS 
PHOTOMETRY OF OUR GC SAMPLE 

2.1 Selection of the GC Sample 

The GC sample u sed in this paper was selected from that of 
Ijiang et al~l 120031) . who published their intermediate-band pho- 
tometry in 13 passbands taken from the BATC Multi-colour Sur- 
vey of the Sky. They estimated the ages using BC96 for 172 M31 
GCs, selected from the Bologna catalogue (Batt istini et al.ll987l) . 
This catalogue, which includes a total of 827 GC candidates, is the 
most comprehensive list of M31 GCs. Of these objects, 353 are 
considered class A (254) and class B (99) objects, i.e., probable 
GCs (or candidates) to a very high and high level of confidence, re- 
spectively; 152 sources are listed as class C, and the others fall into 
the lower-confidence classes D and E. Ijiang et al. I J2OO3I) took all 
candidates in classes A and B as their original GC sample. How- 
ever, only 223 objects overlapped in position with the BATC 58' 
x 58' CCD field 2 . Four o f these are most likely stars rather than 
GCs iBarmbvetaLl fcoOO). An additional 47 GCs were saturated 
in some of the BATC filters, and were therefore not included in the 
final Jiang et al. sample. As a con sequ ence, there are 172 GCs in 
their final sample (see details from Jiang et al. 2003). For the pur- 
pose of est imating accurate GC ages, we selected GCs from the 
Jian g et al .1 <2003l) sample for which the metallicity and reddening 
values had been estimated accurately and homogeneously in previ- 
ous studies. iBarmbv et al. N200ol) and lPerrett et al. 1 1200 2) derived 
metallicities of, res pectively 61 and mor e than 200 GCs and GC 
candidates in M31. lBarmbv etal.N2000l P. Barmby, priv. comm.) 
determined the reddening values for 223 GC candidates. Based on 
cross-id entificati on o f thejjiang et al. |2003) sample and that of 
lBarmbvetai~ll2000l) and lPerrett et al. 1 120021) . we selected 93 ob- 
jects as our initial M31 GC cluster sample. We refer the reader 
to Section 3 for more technical details. When we compared the 
multi-colour photometry of our sample GCs with theoretical stellar 
synthesis models, we could not obtain a reasonable fit to the SED 
of one GC, Bol31. A detailed investigation of the M31 images of 
BATC Multi-colour Survey of the Sky reveals that this GC is su- 
perimposed onto a very high and variable background, so that the 
resulting photometry is highly uncertain. Therefore, we discarded 
this object from our final sample. In ad dition, 037-B327 has been 
studied in detail bv IMa et al. 1 12006b). and is therefore also ex- 
cluded from the present sample. Our final sample in this paper thus 

contains 9 1 GCs. 

Since ljiang et al. Il2003l) 's GC sample selection was only lim- 
ited by the size of their field of view, their GC sample constitutes 
a random selection of M3 1 GCs, provided that there are no signif- 
icant positional biases present in the M31 GC sample as a whole 
between the disc and the outer halo. In order to illustrate this limi- 
tation, Fig. 1 shows the map of all known bona fide M3 1 GCs from 
lGalletietal.l<2004 and the subsample of GCs discussed in this pa- 
per. It is clear that our sample GCs are indeed limited to the central 
region of M31, as determined solely by the BATC observational 
field of view, and without further selection criteria imposed that 
might jeopardise the random nature of our sample selection. How- 
ever, we point out that our sample may be affected by an under- 



1 An SSP is defined as a single generation of coeval stars formed from the 
same progenitor molecular cloud (thus implying a single metallicity), and 
governed by a given initial mass function (IMF). 



2 The BATC field of view covers most of the visible disc of M3 1 ; the con- 
sequence of the limited field of view for the present paper is that we may 
have missed a number of GCs in the extreme outer disc regions along the 
galaxy's major axis, and that the fraction of outer halo GCs is underrepre- 
sented (see aslo Fig. 1). 
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Figure 1. Map of all known bona fide M31 GCs from lGalleti et al. 1 120041) 
and the subsample of GCs discussed in this paper. Open and solid circles 
present the Galleti et al. 1 2004 ) GCs and our subsample, respectively; the 
ellipse traces the D25 isophote of the galaxy. 
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Figure 2. (V — I) distribution of all known bona fide M3 1 GCs from 
Galleti etal.ll2004 and the subsample of GCs discussed in this paper. 



der to show any differences as clearly as possible, we use the 
normalised number of GCs in each sample. Fig. |2| shows that the 
(V — /) colour distribution of our subsam ple of M31 GCs is sim- 
ilar to that of the full Galleti et alj 120041) sample. If anything, the 
disc GCs, by which our subsample is dominated, may be somewhat 
redder than the mean of the entire M31 GC sample. This is not sur- 
prising in view of the possible selection biases discussed above, 
while the main galactic disc of M3 1 may also give rise to redder 
GC colours owing to the presence of fairly significant dust lanes 
(particularly closer to the galactic centre). 

In Section l4~4l we will f urther consid e r to w hich extent our 
GC sample represents that of Ijiang et al. I 120031) . We also note 
that although mo st present ca talogues of M31 GCs are not com- 
plete, [FusTPeccnyt^ G993|) have shown that the M31 GC pop- 
ulation (including the sample in the Bologna catalogue on which 
the present sample is based) is currently fairly complete down to 
V = 18 mag (My ~ —6.5 mag). As such, we believe that we 
have not introduced additional photometric biases into our sample 
selection by using the Bologna catalogue as our master list of M31 
GCs. 



2.2 Intermediate-band Photometry of Our Sample GCs 

The large-field multi-colour observations of M31 were obtained 
in the BATC photometric system, obtained with the 60/90cm f/3 
Schmidt telescope at the Xinglong Station of the National Astro- 
nomical Observatories of China. A Ford Aerospace 2048x2048 
CCD camera, with a 15 /im pixel size is mounted at the Schmidt 
focus. This set-up provides a CCD field of view of 58' x 58', 
with a pixel size of l'.Y. The filters were designed specifically to 
avoid contamination f rom the brig htest and most variable night sky 
emission lines. Ijiang et al. I 120031) extracted 123 images of M31 
from the BATC survey archive, taken in 13 BATC filters at wave- 
lengths from 4,000 to 10,000A between 1995 September and 1999 
December. They combined multiple images through the same fil- 
ter to improve the image quality, and determined the magnitudes 
of 172 GCs in the 13 BATC filters based on the combined images 
using standard aperture photome try, i.e., in esse nce by employing 
the PHOT routine in DAOPHOT IStetson |[l987l) . The BATC pho- 
tometric system calibrates the zero-magnitude level in a similar 
fashion to the spectrophotometric AB magnitude system. For flux 
calibration, the Oke-Gunn primary flux stand ard stars HP 19445, 
HD 84937, BD +26°2606, and BD +1 7°4708 lOke & Gunnll983l) 
were observed on photometric nights lYanetallb oOO). Table 2 of 
Jian g et al. 1 l2003l) lists the resulting BATC photometry of their 172 
GCs, including the 91 GCs in this paper. 



representation of young (< 1 — 2 Gyr-old ) blue luminous compact 
clusters (BLCCs). Fusi Pecci et al"ll2005t) found that such BLCCs 
seem to avoid the inner regions of the M3 1 disc, and appear to be 
clearly projected onto the spiral-arm structure in the outer disc (pre- 
dominantly along the major-axis direction; see their Fig. 5), with 
which they also share kinematic properties (see also IPuzia et al. I 
2005). In addition, there is evidence for a roughly 25 to 50 per 
cent solar- metallicity ([Fe/H] = —0.6 to —0.2 dex ) population of 
older GCs in the M31 halo (e.g. lPuzia et al. 1l200l and references 
therein), while the more metal-poor clusters appear to be sca ttered 
throughout the entire disc-halo syste m (e.g. iBarmbv et al. lEoOOt 
IPerrett et al. l2002UPuzia et al. 120051) . 

F ig. |2] compares the (V — I) colour distributions of the 
iGalleti et aU 120041) M31 GCs with that of our subsample. In or- 



2.3 Broad-band Photometry of Our Sample GCs 

In order to estimate the ages of our sample GCs accurately, we try 
to use as many photometric data points covering as large a wave- 
length range as possible (cf. Ide Griis et al. lEo03bt lAnders et al. I 
2004). Using the 4-Shooter CCD mosaic camera and the SAO in- 
frared image r on the 1.2m te l escop e at the Fred Lawrence Whipple 
Observatory, IBarmbv et al. I l2000l) presented optic al and infrared 
photo metric data for 285 M31 GCs (see Table 3 of IBarmbv et al. I 
2000). For most of our sample GCs, photometri c measureme nts in 
the UBVRI bands were published by Barmbv et al. (2000). For 
the remaining GCs, we refer to Gall eti et al. I 12004I) . who updated 
the Bologna Catalogue with the homogenised optical (UBVRI) 
photometry collected from the most recent photometric references 
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available in the lit erature. In this revised Bologna Catalogue of M3 1 
globular clusters, iGalleti et al. I 120041) did not include the photo- 
metric uncertainties. Therefore, we refer to the original works as 
indicated in Gal leti et al. I 12004 ). For a number of GCs, for which 
U -band photometric uncertainties were not presented in the origi- 
nal references, a photometric uncertainty of .08 mag was adopted 
(see lGalleti et al. 12004) . IGalleti et al. 1 1 2004) took the photometry 
in the UBVRI bands from Barmbv et al. (2000) as a reference 
to obtain a Master Catalogue of photometric measurements com- 
piled in as homogeneously a fashion as possible. All other cata- 
logues containing UBVRI photometry were transformed to this 
reference by applying the offsets derived from objects in common 
between the relevant catalogue and the list o flBarmbv et al. I I2OOO). 
Therefore, the measurements are internally consistent. 



2.4 2MASS Photometry of Our Sample GCs 

As pointed out bv lWorthevNl994l) . the age-metallicity degeneracy 
in optical broad-band colours is Aage/AZ ~ 3/2. This implies 
that the composite spectrum of an old stellar population is indis- 
tinguishable from that of a younger but more metal-rich population 
(an d vice versa) (see also iMacArthur et al. 120041) . Ide Jong 
and I Anders et al. I 120041) showed that this degeneracy can be par- 
tially broken by adding infrared photometry t o the optical colours , 
depending on the age of the stellar population. lCardiel et al. 1 120031) 
found that the inclusion of an infrared band can improve the pre- 
dictive power of the stellar population diagnosti cs by ~ 3 times 
compared to using optical photometry alone. IWu et al.1 120051) 
also showed that the use of near-infrared colours can greatly con- 
tribute to break the age-metallicity degeneracy. In this paper we 
add 2MASS photometry to our broad and intermediate-band op- 
tical photometry to estimate the GC ages to the highest possible 
accuracy. 

Using the 2MASS database, IGalleti et all 120041) identified 
693 known and can didate GCs in M31, and listed their 2MASS 
JHKs magnitudes. IGalleti et all 120041) trans formed all 2MASS 
magn itudes to the CIT photometric s ystem {E lias et al"l 1 1 9821 
1983) using the colour transformations in Carpenter (2001). How- 
ever, we need the original 2MASS JHKs magnitudes for our 
sample GCs in order to compare our observational SED to the 
SSP models, so we reversed this tran sformation using the same 
procedures. Since Galleti et al. (2004) did not give the photomet- 
ric 2MASS JHKs uncertainties, we obtained photometric uncer- 
tainties by comparing the photometric magnitudes with Fig. 2 of 
ICarpenter et al. ifcpOlh . who plot the observed photometric rms un- 
certainties in the time series as a function of magnitude for stars 
brighter than their observational completeness limits. In fact, the 
photometric uncertainties adopted do not affect our results signifi- 
cantly, as we will show below (see section 4.3 for details). 

Finally, we point out that because of the very long wavelength 
coverage from U to the near-infrared, and the large number of pho- 
tometric data points for each GC, our age estimates, and in par- 
ticular the associated uncertainties are highly robust. In fact, de- 
pending on the (internal) extinction curve adopted, they are com- 
petitive with respect to spectroscopic age determinations based on 
either a short wavelength range or a limited number of age diagnos- 
tics (where small-scale differences in the continuum level may have 
significant consequen ces for t he age determinati ons), as d iscussed 
extensively inlde Griis et al.l l2003bl) and Ide Griis et al. I J2005a). 
ISchweizer et al. I l2^04jr provide convincing support for this state- 
ment by showing that spectroscopic age determinations are not nec- 



essarily better or more accurate than ages obtained from photome- 
try. 



3 METALLICITIES AND REDDENING 

To estimate the ages of our sample GCs accurately, we required that 
our GC sample have both independently determined spectroscopic 
metallicities, and reddening values. We used two homogeneous ref- 
erence sources of spectr oscopic metallicities. lBarmbv et al. 1 12001) 
and lPerrett et al. 1120021). 

Barmb v et al. I yOOO) present metallicities of 61 GC candi- 
dates, using the Keck LRIS and the MMT Blue Channel spec- 



trographs. With the Keck LRIS, they used a 600 I mm grating 
with a 1.2 A pixel -1 dispersion from 3670-6200 A, and a reso- 
lution of 4-5 A. With the MMT Blue Channel, they used a 300 t 
mm -1 grating with a 3.2 A pixel" 1 dispers i on fro m 3400-7200 
A, and a resolution of 9-1 1 A. Barmb^yilJj200C ) comp uted the 
absorption-line indices followin^ lBrodi^&*Huchra Nl99d) : subse- 
quently, they used the Brodie & Huchra 1 1990) spectral index - 
metallicity calibration. Brodie & Huchra ( 1990)'s metallicity cali- 
bration is based on six absorption-line indices measured from inte- 
grated cluster spectra, which provides an [Fe/H] estimate accurate 
to abo ut 15 per cent. 

IPerrett et ai~l 120021) list metallicities for more than 200 GCs 
in M31 using the Wide Field Fibre Optic Spectrograph at the 4.2 m 
William Herschel Telescope, which provides a total spectral cover- 
age of ~ 3700-5600 A with two gratings. One grating (H2400B, 
2400 I mm -1 ) yields a dispersion of 0.8 A pixel" 1 and a spectral 
resolution of 2.5 A over the range 3700-^1500 A, and the other grat- 
ing (R1200R, 1200 I mm -1 ) is characterised by a dispersion of 1.5 
A pixel -1 and a spectral res olution of 5.1 A over the range 4400- 
5600 A. |Perrettet aLl<2002l) then calculated 12 absorption-line in- 
dices following lBrocfe & Huchra 1 ll990h . By comparison of the 
line ind ices with published M3 1 GC [Fe/H ] values from pre vious 
studie s iBonoli et al. ll987l : lBrodie & Huchra 1 199rAlBarmbv et al. I 
2000 ), they applied a linear least-squares minimisation. Final clus- 
ter metallicities were determined from an unweighted mean of the 
[Fe/H] values calculated from the CH (G band), Mgfe, and Fe53 line 
strengths. 

We draw the reader's attention to the fact that the M31 GC 
metallicities thus o btaine d were based on the [Fe/H] calibration of 
iBrodie & Huchra lll99Cj) . who included only old GCs in their cali- 
bration sample. Since in this paper we estimate the ages of our M3 1 
GC sample, and find that some of them are young or of intermediate 
age, our age determination may be somewhat biased. However, we 
point out that (i) the main features in the M3 1 GC age distribution 
that we will derive in Section(s) 4.5 (an d 4.6) are similar to those 
derived by authors who do not rely on the[B rodie & Huchra 11 1990l) 
metallicity calibration (e.g. Beas lev et al. 120051) - on a one-to-one 
basis, the ages of most of the GCs we have in common with these 
authors do not differ by more than 1 — 2a, and (ii) we do not have 
the spectroscopic data required to quantify this effect in detail; in- 
stead, we refer to lBeaslev et al. I l2005t) for an in-depth discussion 
on this issue. 

For the reddeni ng values of the sample GCs we refer to 
iBarmbv et al. I l200fj) . who determined the reddening for each in- 
dividual cluster using correlations between optical and infrared 
colours and metallicity and by defining various "reddening-free" 
parameters using their large database of multi-colour photometry. 
IBarmbv etai~l |2p00) found that the M31 and Galactic GC extinc- 
tion laws, and the M31 and Galactic GC colour-metallicity rela- 
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tions are similar to each other. They then estimated the redden- 
ing to M31 objects with spectroscopic data using the relation be- 
tween intrinsic optical colours and metallicity for Galactic clus- 
ters. For objects without spectroscopic data, they used the relation- 
ships between the reddening-free paramet ers and certain in trinsic 
colours based on the Galactic GC data. lBarmbv et al, I fcOOOh com- 
pared their results with those in the literature and confirmed that 
their estimated reddening values are reasonable, and quantitatively 
consistent with previous determinations for GCs across the entire 
M31 disc. In particular. lBarmbv et al. I i200(f) showed that the dis- 
tribution of reddening values as a function of position appears rea- 
sonable in that the objects with the smallest reddening are spread 
across the disk and halo, while the objects with the largest redden- 
ing are concentrated in the galacti c disk. The refore, we adopted the 
reddening values from lBarmbv et all (2000. also P. Barmby, priv. 
comm.) for our GC sample. 



4 AGE DETERMINATIONS OF THE SAMPLE GCS 

4.1 Stellar Populations and Synthetic Photometry 

In evolutionary synthesis models, SSPs are modeled by a collec- 
tion of evolutionary tracks of stars with different initial masses and 
chemical compositions, and a set of stellar spectra at different evo- 
lutionary stages. To estimate the ages of our sample GCs we com- 
pare their SEDs with the updated SSP models of BC03. BC03 pro- 
vide the evolution of the spectra and photometric properties of SSPs 
for a wide range of stellar metallicities. The model set includes 26 
SSP models (both of high and low resolution) based on the 1994 
Padova evolutionary tracks, 13 of which were computed using the 
IChabrierl {2003) IMF assuming lower and upper mass cut-offs of 
mL = 0.1 Mq and mo = 1 00 Mr, r espectively, while the other 
13 were computed using the Salpeter 1 1955) IMF with the same 
mass cut-offs. In addition, BC03 provide 26 SSP models using the 
2000 Padova evolutionary tracks. However, as they point out, the 
2000 Padova models tend to produce worse agreement with ob- 
served galaxy colours. These SSP models contain 221 spectra de- 
scribing the spectral evolution of SSPs from 1 x 10 yr to 20 Gyr. 
The evolving spectra include the contribution of the stellar compo- 
nent at wavelengths from 91 A to 160/iin. In this paper, we adopt the 
high-resolution SSP models co mpute d using the 1994 Padova evo- 
lutionary tracks and a lSalpeterl fl 9551 IMF 3 . We note that although 
the current best estimate of the age of the Universe is of order 13.7 
Gyr, the SSP models and the stellar evolutionary tracks that form 
their basis have been calculated for ages up to 20 Gyr. It is not 
straightforward to correct for this discrepancy; one would need to 
recalculate all stellar evolutionary tracks for all metallicities. This 
situation is exacerbated by the fact that (i) young open clusters, up 
to ~ 5 Gyr, have been used to properly and robustly constrain the 
stellar evolution at younger ages, and (ii) at older ages the obser- 
vational diagnostics are less sensitive to changes in age, resulting 
in significantly greater uncertainties at those ages. These are there- 
fore issues that one needs to keep in mind in the context of physical 
parameters derived for the older GCs. 

Since our observational data are integrated luminosities 

3 We note that because of the slow SED evolution of SSPs at ages in excess 
of a few Gyr, all of the most commonly used spectral synthesis models agree 
very well at these ages. Therefore, the choice of IMF is only important 
for estimating the photometric mass of the cluster, and does not affect the 
determination of the ages of old GCs. 



through our set of filters, we convolved the BC03 SSP SEDs with 
the BATC intermediate-band, broad-band UBVRI and 2MASS 
filter response curves to obtain synthetic optical and near-infrared 
photometry for comparison. The synthetic i filter magnitude can 
be computed as 

rrii = -2.5 log ^ ^^-48.60 , (1) 

where Fx is the theoretical SED and tpi the response curve of the 
i th filter of the BATC, UBVRI and 2MASS photometric systems. 
Here, F\ varies with age and metallicity. 

4.2 Fit Results 

We use a \ 2 minimisation test to examine which BC03 SSP models 
are most compatible with the observed SED, following 

X =2^ -2 . 

where m™° d (t) is the integrated magnitude in the i th filter of a the- 
oretical SSP at age t, m^ bs presents the observed integrated mag- 
nitude in the same filter, and 

^ % ^obs,z "f" ^mod.j ■ (-') 

Here, a^ hs ^ is the observational uncertainty, and <J^ nod t is the 
uncertainty associat ed with the model itself, for the i th filter. 
ICharlotetal.U l996) estimated the uncertainty associated with the 
term crj^ od 4 by comparing the colours obtained from dif ferent stel- 
lar evo luti onary tracks and s pectral libraries. Following W u et al. I 
<2005h and lMa et al.H 2006b1. we adopt a^ od i = 0.05. 

The BC03 SSP models include six initial metallicities, 0.0001, 
0.0004, 0.004, 0.008, 0.02 (= solar), and 0.05. Spectra for other 
metallicities can be obtained by linear interpolation of the appro- 
priate spectra for any of these metallicities. For our sample GCs, 
whose metallicity and reddening values were published by other 
authors, the cluster age is the sole parameter to be estimated (for 
a given IMF and extinction law, which we assume to be universal 
among our GC sample and for the main calibrators). The values 
for the extinction coefficient, R\, are o btained by interpolating the 
interstellar extinction curve of ICardelli etail Jl989h . In Fig.[0 wc 
show the observational SEDs of a representative, randomly selected 
subsample of M31 GCs, and the SEDs of the best-fitting models. 
The fit results for our entire sample are listed in Tables 1 and 2. 

4.3 Check of the importance of the photometric uncertainties 

We will explore and quantify to what extent the photometric un- 
certainties affect our results. To this end, we randomly selected a 
number of GCs, while making sure that they spanned our entire 
age range. The results for a small subsample of our GCs are listed 
in Table 3. From the left to the right column on the same line, the 
photometric uncertainties vary from their original values (1), half 
the original (2), and twice the original (3). From this initial com- 
parison it is clear that the photometric uncertainties adopted do not 
affect our results significantly, in view of the large uncertainties in 
the derived ages, but we note small differences in the best-fitting 
ages towards the high-age end. In Fig.|4]we show the results for the 
entire sample. The subscripts, 1, 2 and 3, refer to conditions (1), 
(2) and (3) above, respectively. A comparison of Figs.|4^ and c on 
the one hand, and b and d on the other show that (i) the best-fitting 
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Figure 3. Observational SEDs of a representative, randomly selected subsample of M3 1 GCs, and the SEDs of the best-fitting models. 



age estimates are essentially independent of the photometric uncer- 
tainty assumed for ages up to about 8 Gyr; (ii) for older ages, the 
best-fitting ages tend to be overestimated by < 1.5 Gyr, in the ma- 
jority of cases, if the photometric uncertainties are halved, and (iii) 
underestimated by < 2 Gyr, for most objects, if the photometric un- 
certainties are doubled. We also note, however, that the actual un- 
certainties on the derived ages are well in excess of these small sys- 
tematic effects at old age; for this reason (and for improved clarity) 
we have omitted the error bars in Figs.|4j; and d. Nevertheless, we 
will need to keep these small systematic offsets in mind when as- 
sessing the absolute uncertainties in our best-fitting age estimates, 
which we will do in the next section. 

4.4 Comparison of age determinations 

Ijiang etal.li2003l) estimated the ages of their 172 M31 GCs, based 
on only the BATC data and on the SSP mode ls of BC96; the 
reddening values ado pted by Ijian^^] |2003j) were also from 
iBarmbv et aTl <200fl . However. Ijiang et al. I <2003l) only used the 
BC96 models for three metallicities, i.e., 0.0004, 0.004 and 0.02, 
and did not linearly interpolate to find the best-fitting metallicities. 
For old GCs, the age/metallicity degeneracy becomes important, so 
it is reasonable to postulate that Jiang et al. N2003I) did not estimate 
the ages of their GC sample as accurately as we have done in this 
paper. Figure[5]shows a comparis on of the age dete rmination for the 
M31 GCs in common between Ijiang et al. U2003I) and this paper. 
It is immediately clear that there is a significant difference between 
I Jiang et all 1 2003) and our results. 

In order to have confidence in our age determinations in this 




5 10 15 20 5 10 15 20 



Age, (Gyr) Age, (Gyr) 

Figure 4. Comparison of the best-fitting GC ages and their dependence on 
the photometric uncertainties assumed. Subscripts 1, 2 and 3 refer to the 
assumption of, respectively, the actual photometric uncertainties, half and 
twice these uncertainties. We have omitted the error bars in panels c and d 
for reasons of clarity; the dotted lines indicate the loci of equality. 



Table 1. Fit Results for our 91 M31 GCs. 



New age estimates ofM31 globular clusters 7 

Table 2. continued 



ID 


Age 


Metallicity 


2 

-^min 




(Gyr) 


[Fe/H] 


(per degree of freedom) 


Bo009 


16.75±4.90 


-1.57±0.26 


2.87 


Bo015 


1.43±0.04 


-0.35±0.96 


1.86 


Bo017 


2.30±0.07 


-0.42±0.45 


0.59 


Bo023 


20.00±4.94 


-0.92±0.10 


2.69 


Bo025 


8.25±2.08 


-I.46i0.13 


1.62 


Bo027 


17.50±2.95 


-1.64±0.32 


3.80 


Bo029 


14.25±4.25 


-0.32±0.14 


5.51 


Bo030 


5.25±1.30 


-0.39±0.36 


3.81 


Bo031 


2.60±0.92 


-1.22±0.40 


4.56 


Bo034 


20.00±4.69 


-1.01±0.22 


3.46 


Bo036 


20.00±4.34 


-0.99±0.25 


2.13 


Bo038 


6.00±0.82 


-1.66±0.44 


0.43 


Bo039 


4.75±0.47 


-0.70±0.32 


0.61 


Bo041 


19.50±5.03 


-1.22±0.23 


0.76 


Bo042 


2.60±0.14 


-0.78±0.31 


2.11 


Bo044 


8.25±2.26 


-1.14±0.37 


1.51 


Bo051 


11.25±2.51 


-1.00±0.13 


1.14 


Bo054 


1.28±0.03 


-0.45±0.17 


3.00 


Bo056 


3.50±0.78 


-0.06±0.10 


3.54 


Bo057 


6.75±1.95 


-2.12±0.32 


1.23 


Bo059 


20.00±1.01 


-1.36±0.52 


1.83 


B0O6I 


7.25±0.87 


-0.73±0.28 


1.73 


Bo063 


9.75±2.58 


-0.87±0.33 


1.43 


Bo064 


18.00±2.47 


-1.55±0.30 


1.85 


B0O68 


1.80±0.04 


-0.29±0.59 


1.33 


Bo073 


7.25±0.87 


-0.64±0.46 


1.78 


Bo076 


1.68±0.10 


-0.72±0.06 


1.78 


Bo082 


6.00±0.87 


-0.80±0.18 


3.35 


B0O86 


15.25±2.72 


-1.74±0.17 


1.25 


B0O88 


7.75±0.84 


-I.81i0.06 


0.60 


Bo093 


20.00±5.44 


-I.03i0.12 


1.64 


Bo094 


3.25±0.30 


-0.17±0.45 


1.34 


Bo096 


9.00±1.65 


-0.26±0.43 


1.81 


Bo097 


2.30±0.23 


-1.21±0.13 


1.19 


Bo098 


3.00±0.11 


-0.67±0.58 


1.41 


Bo 102 


1.80±0.24 


-1.57±0.10 


2.05 


Bo 1 03 


7.75±1.15 


-0.56±0.62 


1.17 


Bo 1 05 


10.25±1.82 


-I.13i0.32 


0.56 


Bo 106 


20.00±7.11 


-0.86±0.68 


2.18 


Bo 1 07 


5.25±1.84 


-1.18±0.30 


2.25 


BollO 


8.00±2.43 


-1.06±0.12 


2.00 


B0II6 


7.50±0.76 


-0.88±0.12 


0.67 


Bol25 


4.00±0.71 


-1.52±0.08 


1.62 


Bo 127 


9.75±3.03 


-0.80±0.14 


3.88 


Bol30 


6.50±1.74 


-1.28±0.19 


1.71 


BoI34 


5.50±2.28 


-0.64±0.08 


5.47 


Bol35 


20.00±2.10 


-1.62±0.04 


1.18 


Bol37 


5.50±1.08 


-1.21±0.29 


1.39 


Bol43 


4.75±0.76 


0.09±0.42 


1.99 


Bol48 


7.50±2.34 


-1.15±0.34 


3.07 


Bol49 


4.25±1.11 


-1.35±0.25 


3.38 


Bol51 


12.50±2.61 


-0.75±0.18 


2.05 


Bol52 


17.50±4.76 


-0.87±0.49 


2.12 


Bol53 


8.00±0.83 


-0.08±0.33 


1.10 


Bol54 


1.61±0.17 


-0.45±0.63 


4.36 


Bol58 


20.00±5.57 


-1.02±0.02 


3.06 


B0I6I 


20.00±1.00 


-1.60±0.48 


3.64 


Bol63 


11.75=1=1.60 


-0.36±0.27 


1.87 


Bo 164 


1.61±0.18 


-0.09±0.40 


1.88 


Bo 165 


11.50±3.08 


-1.80±0.32 


1.49 



ID 


Age 


Metallicity 


y 2 ■ 




(Gyr) 


[Fe/H] 


(per degree of freedom) 


dOID / 


^ ^H-U 1 1 7 

j. jUr_t l . IZ 


n /i7_i_n 77 
-u.4Zr_tu.zj 


1 1 s 

1 . 1 J 


dOI / 1 


1 fl 7^-U 1 Q7 

1U. / Jit 1 .0 / 


n /i 1 _ i— n n/i 

-U.41 _tU.U4 


7 AO 

z.4y 


dOI /4 


1 fi ^n_i_ 1 Qfi 


1 A7_l_n 77 

-i.o/±u.z / 


1 77 
l.Z / 


TJ r\ 1 7 fi 

DOI /o 


1 ^ <n_i_7 77 

1 J. JUltZ. J 1 


i «i _j_n 1 7 
- 1 .3 1 ztU. J Z 


1 D7 
1 .UJ 


T3r. 1 70 

001 /y 


O 7^_!_7 1 7 

y. / j±z. 1 / 


-1 . iUztU.UZ 


n 07 

u.y / 


DOloU 


fx 7^-1- 1 77 
O.Zjlt 1 . J / 


-i . lyrxu.u / 


1 .0 / 


0OI0Z 


1 A 7 ^ -U 7 7 

lo. / JirZ.yj 


1 HA _!_n 1 7 
-l.Z4±U. IZ 


1 . J J 


T3 r\ 1 fi 7 


7 nn-tn no 
j.uuztu.uy 


n i o_i_n q 1 
-u. iy_tu.ji 


1 .UJ 




1 AQ_l_n no 
i.OoirU.Uo 


n i7_i_n /in 
-U.J / ±U.4U 


1 £7 
l.DJ 


T3<-.1 fi^ 


Q 7<_l_7 7/1 
y. / J_T_:J.Z4 


n 7^-un ns 

-U. /DztU.Uo 


7 £7 
Z.O / 


d „ 1 on 


1 7 7<_l_7 7,1 
IZ. / J±Z.Z4 


i nu_n no 
- 1 .ujztu.uy 


n ^7 
U.j / 


T3n 1 Q7 

ooiyj 


1 ^ ^n~l-7 07 

1 j. j\j± j.yz 


-U.44ztU. J / 


1 AQ 

i .4y 


T2 1 Q7 


U.ZjrxU.Uj 


"U.4j±U. JO 


7 A7 
Z.O / 


Ti^.7n 1 

DOZU1 


1 n 7<_i_/i 7A 

1U. / J_±4.Z0 


-1 .Uur_tu.zi 


7 A/I 
Z.04 


oOZUj 


O.ZjirU.oO 


n QH-Ln 77 

-u.yuztu.jz 


1 A ^ 

1 .4j 


DOZU4 


7 nfl-Ufl Q7 

/.uuitu.yj 


-U.oUztU. J / 


7 D7 
Z.UZ 


dOZUj 


7n nn I 7 7Q 

zu.uu_rz.zy 


1 7/1 _Lfl 1 7 
- 1 . J4±U. J J 


l.OJ 


dOZUD 


7n nn4-7 <n 

ZU.UUr_tZ.jU 


-i .4Jr_tu. iu 


7 in 
j. ju 


Bo209 


20.00±0.74 


-1.37±0.13 


1.78 


B0211 


19.00±1.68 


-1.67±0.52 


1.03 


Bo213 


7.50±1.61 


-1.02±0.11 


1.89 


Bo214 


4.00±1.50 


-1.00±0.61 


2.63 


Bo217 


10.25±2.15 


-0.93±0.14 


1.05 


B0218 


19.00±2.30 


-1.19±0.07 


2.06 


B0220 


6.00±0.83 


-1.21±0.09 


0.96 


B0221 


16.75±3.57 


-1.29±0.04 


1.47 


B0222 


7.75±1.46 


-0.93±0.95 


5.09 


Bo224 


5.50±1.00 


-1.80±0.05 


0.45 


Bo225 


15.50±4.89 


-0.67±0.12 


2.52 


Bo228 


10.00±1.86 


-0.65±0.66 


0.98 


Bo235 


4.00±0.59 


-0.72±0.26 


1.65 


Bo037 a 


12.4±3.2 


-1.07±0.20 


1.43 



Note: a values taken from lMa etaO feo06b). 



Table 3. Check of the effects of the photometric uncertainties. 



Name 




Best-fitting age (Gyr) 




(original errors) 


(0.5x original) 


(2x original) 


Bo015 


1.43±0.04 


1.43±0.04 


1.43±0.05 


Bo 137 


5.50±1.08 


5. 75 ±1.09 


4.75±1.09 


Bo228 


10.00±1.86 


10.50±1.72 


9.75±2.11 


Bol78 


15.50±2.37 


16.75±2.36 


14.50±2.32 


Bol35 


20.00±2.10 


20.00±1.34 


18.50il.89 



paper, we must make an effort to understand the main cause of 
these significant differences b etween the age estimates presented 
here and those of ljiang etal.N2003l) . Since we suspected that the 
age-metallicity degener acy may be the pri ncipal culprit in this re- 
gard, we used Table 3 in Jiana et al. 1 2003) to produce Fig.|6| Here, 
we show the difference in age versus the difference in metallic- 
ity for the M31 GCs in common between our sample s; we used a 
straightforward conversion from Z, as listed in ljiang etal.U 2003). 
and [Fe/H], assuming for the sake of simplicity that these parame- 
ters are one and the same. The error bars along the age axis are a 
combination of the uncertainties in age from both th e present pa - 
per and the 20 per cent uncertainties quoted in ljiang et al. I J2003h : 
the error bars along the metallicity axis represent the uncertainties 
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Figure 5. Comparison of age determinations for the M3 1 GCs in common 
between Ijiang et aTl 120031) - based on lBruzual & Chariot! fl99f|) and on 
the BATC multi-band photometry; the error bars correspond to their "typi- 
cal" 20 per cent uncertainties - and this paper. 



in [Fe/H] from Table 2 only. From Fig.[6|it is immediately obvious 
that there is indeed a very strong age-metallicity degeneracy present 
among our GC sample; this is also (and just as significantly so) the 
case if we had used "Age" instead of "log(Age)" as our horizontal 
axis (not shown). The correlation shown in this figure is excellent, 
and in the sense expected if our new age estimates are indeed more 
accurate 4 ; the scatter is li kely partially owin g to the discrete steps 
in metallicity adopted bv ljiang et al. 1 12003|), as well as to the in- 
trinsic uncertainties in the metallicity determination and possibly 
to a residual (but small; see footnote 4) age-metallicity degeneracy 
in the age determinations presented here. In any case, the scatter in 
the relationship is well inside the ~ la level. Thus, here we have 
shown that we fully understand the significant differences between 
our age determinations and those o f ljiangetal.i l 2003). We believe 
to have presented significant improvements in the present paper, as 
well as clear evidence that the previous age determinations of statis- 
tically large samples of M31 GCs were more significantly affected 
by the age-metallicity degeneracy affecting old GCs than our new 
age estimates. 

We will now explore whether our GC sample constitutes a ran- 
dom subsample of the Jiang et al. (2003) GC sample. Fig.0shows 
the age distribution for the full Jiang et al. 1 120031) GC sample, as 
well as for the full GC sample selected in this paper. Both distri- 
butions are based on the age determinations of ljiang et al. I 120031) . 
In Fig.Q the numbers have been normalized, i.e., in each bin the 



4 We note that the tight correlation is predominantly driven by the age- 
metallicity degeneracy evident in the Jiang et al. 1 200 results: the distri- 
bution of data points in the A log(Age) versus [Fe/H] <Jiangetal.l2003T) 
plane shows a clear gen eral trend of dec reasing metallicity with increas- 
ingly large ages in the |jiangetal"1<2003l) results; such a clear correlation 
is not evident in the Alog(Age) versus [Fe/H] plane denned by our new 
data. 
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Figure 6. Age and metallicity differences between the present paper and the 
determinations by Jiang et al. • 201).';): see text for details. 



yyyyyy//yy\ GC subsample of this paper 

| | Full GC sample of Jiang et al. (2003) 




Figure 7. Age distributions of the M31 GCs in ljiang et al. I l2003h a nd in 
this paper, both based on the ages determined by Jiang et al. 1 2003J). The 
age distribution of the GCs discussed in this p aper can therefore immedi- 
ately be compared to the full distribution of the lJiang et al. 1 12003) sample. 



number of objects is divided by the total number of objects. To first 
order, Fig.|7]sho ws that our sample c onstitutes a randomly selected 
subsample of the lJiang et al. 1 120031) GC sample; the main features 
of the age distribution are apparent in both distributions. Any ap- 
parent differences, such as the relative ratio of the two age peaks 
at ages below 10 Gyr, can be attributed to Poissonian uncertain- 
ties at the ~ 1 — 2a level. To quantify the similarities and differ- 
ences between both age distributions, we employed a Kolmogorov- 
Smirnov (KS) test. The maximum value of the absolute difference 
is jD max = 0.1259 for these two samples, containing 91 and 172 
data points, respectively. The probability of obtaining this value for 
-Dmax is 28.3 per cent, which means that there is a small, 28.3 per 
cent probability that the two distributions are different. 
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Figure 8. Age distribution of our 91 sample GCs. Poissonian statistical un- 
certainties are shown. 



Figure 9. Relative uncertainties in the M31 GC ages determined in this 
paper, as a function of age. 
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4.5 The M31 GC age distribution 

Fig. [8] shows the (newly determined) age distribution of our 
sample GCs. It shows that there exist young and intermediate- 
age subpopulations among the M31 GC population, peaking at 
ages of ~ 3 and ~ 8 Gyr, respectively. The "usual" com- 
plement of GCs of simil ar age as the o l d Galactic GCs , previ- 
ously discussed bv [ Barmbv et al. | |200fj" 
Burs tein et al. I 120041) . and IPuzia et aUTYlUUa) . is 
present as well. Assuming Poissonian statistical errors, the two 
(young and intermediate-age) subpopulations are indeed significant 
at the ~ 3<r level. In order to assess the robustness of these features 
in the newly derived age distribution of our M31 GC sample, we 
need to take into account the corresponding uncertainties in age, as 
well as the bin size (in Gyr) used to construct the age distribution. 

A comparison of the uncertainties in our derived GC ages 
(obtained from our SED fitting routine) as a function of age (not 
shown), reveals that the errors increase significantly with increas- 
ing age. This is simply owing to the reduced age resolution of older 
broad-band SEDs, since SSP colours are not as sensitive to age vari- 
ations at old ages (i.e., in excess of a few Gyr) compared to their 
earlier evolution. For GC ages up to ~ 5 Gyr, the uncertainties in 
our derived cluster ages are generally < 1 Gyr, in most cases sig- 
nificantly so. Both the spread and the upper limit of the age uncer- 
tainties increase, the upper limit roughly linearly. The latter reaches 
~ 1.5 Gyr at an age of ~ 7 — 8 Gyr. This is, in fact, the age range of 
interest for our assessment of the robustness of the doubly-peaked 
M31 GC age distribution. Fortunately, as we discussed in Section 
14.31 the derived ages are insensitive to the photometric uncertainties 
assumed at these young and intermediate ages. We therefore adopt 
the uncertainties on the derived ages as the basis for our bin sizes; 
the doubly-peaked age distribution remains a significant feature at 
the ~ 3<t level for age bins of up to 1.5 Gyr. We conclude, there- 
fore, that this feature is real and robustly detected. In addition, in 
Fig. 9 we plot the relative uncertainties in our derived GC ages as 
a function of age. It is clear that, although the absolute age errors 
increase significantly with increasing age, the relative uncertainties 
are rather similar over the full range of derived ages, with a mean 
of AAge/Age ~ 0.17 ± 0.10 (Iff spread). 

The age peak at ~ 7 — 8 Gyr appears to roughly coincide 
with the intermediate-age GC population aged between ~ 5 and 8 



Gyr discussed bv lPuzia et al. I l2005h . which they found to have a 
mean metallicity of [Z/H] ~ —0.6 dex (although they also note that 
their sample is biased towards high-metallicity objects). It is worth 
noting that a similar intermediate-age star-forming event may have 
affected the resolved stel lar population, at least along M31's minor 
axis: Brown et al .1 120031 found indications of an intermediate-age 
resolved stellar population, of ~ 6 — 8 Gyr and with [Fe/H] > — 0.5 
dex, from their HSTI ACS observations. 

The young-age peak at ~ 3 Gyr we detect in this paper has not 
been discussed before; in view of the small age uncertainties (see 
above) at young ages, it is unlikely that this populatio n represent s 
the < 2 Gyr-old population of BLCCs of lFusi Pecci et al. I J20051) . 
Instead, we argue that there may have been an additional violent 
star-forming event that triggered the formation of a GC subpopula- 
tion in the disc of M3 1 some 3 Gyr ago. 

It is well-known that the abundance distributions of GCs in 
many gala xies are bimodal, including that of the M31 GC sys- 
tem (e .g. iBarmbv et al. 1120 00: Perret t et al. 1120021) . iForbes et al. I 
1 1997) found that the metal-rich GCs in elliptical and cD galax- 
ies are closely coupled to the resolved stellar populations of their 
parent galaxies, but the metal-poor GCs are largely independent of 
the galaxies. They concluded that the metal-poor GCs were formed 
at the time of galaxy formation, and the metal-rich GCs at a later 
stage. If so, this may imply that the age distribution of GCs should 
be bi- or multi-modal. Despite our incomplete coverage of the M31 
halo and the extreme outer disc regions along the galaxy's ma- 
jor axis, based on Fig. [5] we can confirm that, in general, the age 
distribution of the M31 GCs is u nlikely to be best described by a 
monomodal distribution (see also lPuzia et al. lEo 05. for an exten- 
sive discussion, based on their high-metallicity biased GC sample). 

The age distributions of both the full M3 1 GC sample and of 
the field stars, combined with the detecti on of a metal-poor GC 
population obeying thin-disc kinematics ( Mor rison et al. 1 12004) . 
poses serious problems for our understanding of the formation and 
evolution of the galaxy's disc. On the one hand, the tightly con- 
strained kinematics argue for a relatively quiescent thin-disc evo- 
lution since early times (but see lAbadi et al. lEo03l) . while on the 
other hand the formation of massive star clusters requires violent 
con ditions (such as galaxy mergers) ra ther than in situ formation 
(see lFusi Pecci et al. l2005l : IPuzia et al. I2005L for a discussion). 
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Table 4. Mean age-metallicity relation for the M31 GCs. 
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Figure 10. Distribution of the ages and metallicities of our 91 M31 GCs. 
EiTor bars have been omitted for reasons of clarity. The filled circles with 
their error bars represent the mean distribution and the Gaussian c's of the 
data points; the age ranges used to obtain the mean values are indicated at 
the bottom of the figure by horizontal bars. 



4.6 A Relationship between Age and Metallicity? 

The Galactic GC system is well-known to obey a clear 
age-metallicity relationship, in the sense of older GCs be- 
ing c haracteris ed by generally lower metallic it ies (e.g. 
ISaraiedini & Christopher]! 19891 IChabover et al. Hl99fJ) . iMa et al. I 
(2006b) simultaneously obtained ages and metallicities for 33 
M31 GCs and GC candidates, based on the BC03 SSP models, 
but did not find any significant relationship between their ages and 
metallicities. However, because o f the age-metallicit y degeneracy 
affecting the analysis done to the Ma et al. (2006b) GC sample, 
the ages and metallicities could not be obtained sufficiently 
accurately. Therefore, using the ages and metallicities of the M31 
GCs discussed in this paper, we re-analysed the M31 GC sample. 
Fig. Ei| shows the metallicity as a function of age for the GC 
sample discussed in this paper. It appears that a general trend 
between GC age and metallicity may exist, in the sense that the 
metal-poorer clusters are older, although the scatter in metallicity 
at a given age is significant. It may be better to postulate instead 
that there is a general lack of both low-metallicity younger GCs 
and of the highest-metallicity clusters at the oldest ages. Error bars 
have been omitted for reasons of clarity. 

In order to assess the robustness of these statements, we pro- 
duce the run of the mean of the data points across the age range, 
shown by the filled circles with the error bars. The error bars in- 
dicate the Gaussian er's of the data points in each age range, indi- 
cated by the horizontal bars at the bottom of the figure. Table 4 is a 
compilation of the data points obtained for the mean values shown 
in Fig. 1101 The trend outlined by these mean data points confirms 
the statements made here regarding the possible existence of an 
age-metallicity relation among the M31 GC sample. In view of our 
carefully validated random sampling of our M3 1 GCs, we believe 
that this age-metallicity relation is representative of the entire M31 
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GC population, providing that our positional (field of view) restric- 
tions did not invalidate the random sampling (see Section 2.1). In 
view of the latter concern, we point out that our GC sample may 
be affected by an underrepresentation of BLCCs; these objects are 
mostly aged between ~ 500 Myr and ~ 2 Gyr, and ma y in fact have 
metallicities in the range -1.0 < [Fe/H] < 0.0 <Fusi Pecci et al. I 
2005). If BLCCs are indeed underrepresented in our sample, their 
additional inclusion would, in fact, strengthen the age-metallicity 
relationship revealed in this paper. 

A similar, al though somewhat less robust result was found by 
Jiang et al. 1 <2003l) . who reported to have determined that the more 
metal-poor M31 GCs appear to be older than their more metal-rich 
co unterparts. This is, in essen ce, a similar conclusion as reached 
by iBarmbv & Huchral l2000l) and iPuzia et al. I J2005I) . while the 
combination of positional information and the age-metallicity rela- 
tionship found here is supportive of the existenc e of a sma l l radia l 
metallicity gradient in M31, discusse d by Bar mbv et al. I (2000); 
IPerrett et all fc002h ; |Puzia et al. 1 120051) and references therein. 

In Section 4.2, we estimated the ages for our 91 GCs us- 
ing a x 2 minimisation method, in which we inclu ded the (system- 
atic) model uncertainties following Ch ariot et al. 1 1 1996). We added 
these in quadrature to the observational uncertainties, in a similar 
fashion as lWuetal.1 (2005) and IMa et al"l <2006bl) . However, we 
should be aware that this procedure could potentially somewhat 
confuse our statistical results, since the model uncertainties might 
introduce biases in the age-metallicity correlation, while the ob- 
servational uncertainties smooth and puff-up such correlations. On 
the other hand, we point out that the model uncertainti es are of the 
same o rder as the measurement uncertainties quoted bv ljiang et al. I 
(2003), at least for a large fraction of their measurements, so that 
the effect of adding the model uncertainties in quadrature implies, 
in essence, less than a doubling of the uncertainties, similar to what 
we did in Section 4.3. As such, we are confident that we have not 
introduced major biases into our fitting routine by taking into ac- 
count the uncertainties inherent in the models (note that we only 
increase the uncertainties following this procedure; we did not add 
a systematic offset to the measurements. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we accurately re-determined the ages of 9 1 M3 1 GCs, 
based on improved data, updated theoretical stellar synthesis mod- 
els and sophisticated fitting methods. In particular, we used photo- 
metric measurements of the 2MASS, which can partially break the 
age-metallicity degeneracy, in combination with optical photome- 
try. We showed robustly that previous age determinations based on 
photometric data were affected significantly by this age-metallicity 
degeneracy. 
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Except for one cluster, the ages of our other sample GCs are 
all older than 1 Gyr. Their age distribution shows populations of 
young and intermediate-age GCs, peaking at ~ 3 Gyr and ~ 8 
Gyr respectively, as well as the "usual" complement of well-known 
old GCs, of similar age as the majority of the Galactic GCs. The 
young-age peak at ~ 3 Gyr we detect in this paper has not been 
discussed before; in view of the small age uncertainties at young 
ages, it is unlikely that this population re presen ts the < 2 Gyr-old 
population of BLCCs of Fu si Pecci et al.1l2005h . Instead, we argue 
that there may have been an additional violent star-forming event 
that triggered the formation of a GC subpopulation in the disc of 
M31 some 3 Gyr ago. 

The distributions of the ages of both the GCs and the field 
stars, combined with the existence of a metal-poor GC population 
in M31 obeying thin-disc kinematics, poses serious problems for 
our understanding of the formation and evolution of the galaxy's 
disc. While the tightly constrained kinematics argue for a relatively 
quiescent thin-disc evolution since early times, the formation of 
massive star clusters requires violent conditions (such as galaxy 
mergers) rather than in situ formation. 

Our results also show that although there is significant scat- 
ter in metallicity at any age, there is a noticeable lack of young 
metal-poor and old metal-rich GCs, which might be indicative of an 
underlying age-metallicity relationship among the M3 1 GC popu- 
lation. 
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